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ABSTRACT 



We discuss the different pliysical processes contributing to tlie infrared continuum 
of AGN, assuming that both photoionization from the active center and shocks ionize 
and heat the gas and dust contained in an ensemble of clouds surrounding the nucleus. 
In our model, radiation transfer of primary and secondary radiation throughout a cloud 
is calculated consistently with coUisional processes due to the shock. We consider 
\ that the observed continuum corresponds to reprocessed radiation from both dust 

and gas in the clouds. CoUisional processes are important in the presence of shocks. 
The grains are sputtered crossing the shock front. The models are constrained by 
CO . sputtering as well as by the far-infrared data. The model is applied to the continuum 

C3 ' of Seyfert galaxies from which best estimate of the nuclear, stellar subtracted, emission 

is available. The results show that radiation-dominated high velocity clouds are more 
^ , numerous in Seyfert 1-1.5 whereas shock-dominated low velocity clouds are dominant 

I ' in Seyfert type 2. This result is in full agreement with the unified model for AGN, 

O \ by which high velocity clouds, placed deeper into the central region and therefore 

reached by a more intense radiation, should play a more significant role in the spectra 
^ ' of broad line objects. We could therefore conclude that in type 2 objects, radiation 

is partly suppressed by a central dusty medium with a high dust-to-gas ratio. Once 
the model approach is tested, a grid of models is used to provide a phenomenological 
k>( analysis of the observed infrared spectral energy distribution. This empirical method is 

. a useful tool to rapidly access the physical conditions of the AGN emitting clouds. For 

' this, analytical forms are derived for the two processes contributing to the infrared 

emission: dust emission and thermal bremsstrahlung produced by the NLR clouds. 
Their relative contribution provides a measurement of the dust-to-gas ratio. 
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1 INTRODUCTION 

The presence of a dusty torus around the central engine was 
proposed by Antonucci & Miller (1985) to explain the broad 
emission-lines observed in the polarized light of Seyfert 2 
galaxies. The torus absorbs a large fraction of the nuclear 
emission and reradiates it in the infrared. This is the basis 
of the so-called AGN unification scheme largely analysed in 
the literature and that provides a geometrical explanation 
for the observed differences between broad line and narrow 
hne objects. 

After the pioneering work by Edelson & Malkan (1986), 
the characteristics of the continuum due to dust heated by 
the optical-UV nuclear radiation has been analyzed by sev- 
eral authors (Barvainis 1992, Pier & Krolik 1992; Laor & 



Draine 1993; Granato & Danese 1994; Granato, Danese & 
Franceschini 1997, Nenkova, Ivezic & Elitzur 2002) . 

The infrared emission is generally associated to a dusty 
torus surrounding the central source of the AGN. Although 
proposed by Pier and Krolik (1992,1993), the calculation of 
the radiative transfer radiation throughout a clumpy torus 
has only been recently performed by Nenkova et al. (2002) . 
The AGN obscuring region is modelled as a toroidal dis- 
tribution of dusty clouds placed at different radii from the 
center, each with same optical depth. The optically thick 
clumps are heated by the nuclear radiation as well as from 
the radiation originated in other clumps. The dust temper- 
ature is much higher in their illuminated side. The authors 
claim that heating by diffuse radiation is highly inefficient 
and can be neglected. Following their results, a clumpy torus 
accounts for both, the often observed broad IR bump ex- 



2 M. Contini, S.M. Viegas, & M.A. Prieto 



tending up to 100 nm and the absence of the 10 /im silicate 
feature in type 1 - on axis view- objects. 

In general, two type of models are used in the litersr 
ture: a compact torus (radial dimension smaller than a few 
pc) with a large optical depth, and a more extended (up to 
hundreds pc) moderately thick disk. A common character- 
istic to all is that dust survival in the vicinity of the AGN is 
set by its evaporation limit. Therefore, dust can reach tem- 
peratures up to 1000-1500 K. Overall, these models provide 
a reasonable fit to the observed spectral energy distribution 
(SED) The torus models usually assume a high equatorial 
opacity leading to flat spectra for type 1 and steep spec- 
tra for type 2 galaxies. However, as pointed out by Alonso- 
Herrero ct al. (2003), the large number of spectral indices 
observed in both types of galaxies, as well as the large num- 
ber of intermediate and type 1 objects with steep SEDs, can 
not be explained by the torus model. 

Other simpler descriptions of the IR SED include e.g., 
that by Blain, Barnard, & Chapman (2003) who describe the 
far-IR to submillimcter SED of dusty galaxies with different 
forms of modified black-body models. 

Although most available models provide in a way or 
other a fair representation of the IR SED, we believe that 
a more sounded physical approach is needed. An important 
problem affecting all these models - those mentioned above, 
as well as in the case of Blain's et al. approach, is that the 
results are independent of the dust-to gas ratio. 

For several years, we are pursuing the development 
of self-consistent models aimed at explaining both the 
emission- line and SED of AGN. These models have proved to 
be successful in a number of representative Scyfcrt galaxies 
(Contini, Prieto, Viegas 1998a,b, Contini, Viegas, & Prieto 
2002). In our modelling approach, the infrared-optical con- 
tinuum is due to dust emission and free-free emission from 
the clouds powered by the nuclear radiation, diffuse radi- 
ation and by shocks. Shocks arc due to the cloud motion 
through a dilute gas (Contini & Viegas- Aldrovandi 1990, 
Viegas & Contini 1994, Contini & Viegas 2000). An impor- 
tant aspect considered in that approach arc the conditions 
of dust survival when effects as sputtering and grain-gas 
collisions are taken on board. Thus, these models offer an 
alternative to the torus model, particularly for the objects 
showing infrared emission as extended as the narrow-line 
region. 

Overall, we consider a model to be well constrained 
when both the complete SED and the emission-line spec- 
trum are simultaneously modeled. When possible, that has 
been our approach, e.g. in the references quoted above. How- 
ever, in this paper, we focus on the optical-IR continuum 
only, since in this wavelength range the best estimate of the 
nuclear flux emission and of shock velocities for a reasonable 
number of Seyfert galaxies of different types (the sample by 
Alonso-Herrero et al. 2003) is available. 

The paper is structured as follows. A short review of 
the processes involving dust is given in Sect. 2. Numerical 
simulations accounting for the presence of both shocks and 
photoionization in the NLR are presented and applied to 
the SEDs of Scyfcrt galaxies, described in Alonso-Herrero 
et al. (2003), in Sect. 3. Concluding remarks appear in Sect. 
4. An empirical method derived from the models, aimed at 
providing a first assessment of the nature of the observed IR 
continuum, is presented in the Appendix. 



2 DUST AND GAS EMISSION 

When the physical conditions of a dust cloud is due to a 
combination of photoionization and shock, the dust heating 
and temperature will depend on the UV radiations as well 
as on the post-shock temperature. 

The maximum grain temperature, which would account 
for the near infrared emission, is limited by sublimation. An 
estimate for some Seyfert 2 galaxies shows that the critical 
distance for grain survival, i.e., the distance where the grain 
temperature is equal to the evaporation temperature, for 
both graphite and silicate, is in the range 0.25 to 50 pc, 
depending on the luminosity of the central source (Vaceli et 
al. 1993). In addition, the value of the temperature reached 
by the grains depends on the presence of small grains and 
on a strong UV radiation from the central source. 

Several lines of evidence arc pointing towards the de- 
struction of very small grains in regions of high UV radia- 
tion intensity (Rowan- Robinson 1992). Destruction of dust 
by those effects are relevant when gas velocities about or 
larger than 300 kms~^are present; still, velocities of ^ 100 
kms~^can have an effect on dust destruction as we shall see 
below. 

In the following, the effect of photoionization and shocks 
on the temperature and survival of dust grains in an emitting 
cloud of the NLR is reviewed. 



2.1 Description of our models: the code SUMA 

Our models result from numerical simulations with the code 
SUMA (Viegas & Contini 1994, Viegas & Contini 1997, and 
references therein) which requires the following input pa- 
rameters: the shock velocity, Vg, the preshock density, no, 
the preshock magnetic field. Bo, an ionizing radiation spec- 
trum, the chemical composition of the gas and a dust-to-gas 
ratio d/g. As ionizing continuum, a power-law character- 
ized by the ionizing flux at the Lyman limit Fh (in photons 
cm~^, s~^, eV"'^) and a spectral index, or a combination of 
spectral indices in different wavelength ranges, is considered. 
A plan-parallel symmetry is adopted with the calculations 
starting at the shocked edge and the cloud divided in slabs. 
The number of slabs is chosen according to the precision 
required by the calculation of the physical conditions. The 
geometrical thickness of the cloud, D, is also an input pa- 
rameter in the case where shock and photoionization act on 
the opposite edges of the cloud. This is the common case for 
clouds flowing outward from the nucleus. 

The models discussed here have been obtained assuming 
a spectral index of 1.5, cosmic chemical abundances for He, 
C, N, O, Ne, Mg, Si, S, Ar, and Fe (Allen 1973), and Bo 
= 10~* gauss. The other input parameters are variable in 
ranges established by the observations as well as by previous 
modelling of the NLR. 

As the calculations account for the sputtering in the dif- 
ferent zone downstream of the shock front (Viegas & Contini 
1994), the distribution of the grain sizes along the cloud is 
automatically derived by SUMA starting from an initial size. 
In AGN, carbon and PAH grain features arc relatively less 
prominent than, for example, in starburst galaxies, there- 
fore only silicate grains are considered (Siebenmorgen et al 
2003). 
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Figure 1. The dust reemission peak relative to bremsstrahlung 
emission for a prototype model M3 (Table 1) 



Self-consistent calculations lead to the continuum and 
emission-line spectra which are compared with the data. 



2.2 Description of our models : the IR emission 

As Nenkova ot al. wo also consider an ensemble of dusty 
clouds placed at different distances from the active center. 
In our model, radiation transfer of primary and secondary 
radiation throughout a cloud is calculated consistently with 
coUisional processes due to the shock. Moreover, we consider 
that the observed continuum corresponds to reprocessed radi- 
ation from both dust and gas in the clouds. For NLR clouds 
moving outwards, in the edge of the cloud facing the pho- 
toionizing source, dust and gas are mainly heated by the 
primary radiation, although the effect of the diffuse radi- 
ation from the ionized gas cannot be neglected. Collision 
effects are less important because they increase with tem- 
perature, which is usually smaller in the photoionized region. 
The temperature of the gas, in fact, does not exceed 2-3 10* 
K in this region. 

On the other hand, at the opposite edge, the gas is 
coUisionally heated to relatively high temperatures by the 
shock. The post-shock temperatures depend on the shock 
velocity (oc Vs'^). In these models, dust and gas are coupled 
throughout the shock front. Collision phenomena dominate 
at high temperature leading to relatively high temperatures 
of gas and grains. Diffuse radiation from the hot gas bridges 
between the shock dominated and the radiation dominated 
sides of the cloud. 

The structure of the gas downstream is determined by 
the shock. In fact, the density distribution results from com- 
pression and the temperature distribution is strongly af- 
fected by the cooling rate. Due to gas and dust mutually 
heating and cooling, the maximum temperature of dust de- 
pends on the shock velocity, which therefore determines the 
frequency of the reradiation peak in the infrared. An an- 
alytical relation between the wavelength of the maximum 



intensity of dust reradiation in the IR spectrum and shock 
velocity is given by Draine (1981). 

The hot NLR clouds produce bremsstrahlung. Depend- 
ing on the temperature distribution in the cloud, this mech- 
anism may contribute to the infrared continuum (see, for 
instance, Contini & Viegas-Aldrovandi 1990, fig 2c). The 
intensity of dust IR emission relative to bremsstrahlung de- 
pends on the dust-to-gas ratio d/g. Usually the dust con- 
tribution dominates in the far infrared range but in the 
near-IR, both processes may contribute to the contirmum as 
shown in Fig. 1, where the single-cloud continuum of model 
M3 (Table 1) is plotted. 

The behaviour at high frequency of the free-free emis- 
sion is largely dependent on the gas velocities. In Fig. 1, it 
shows a linear behavior at low frequencies, then a local max- 
imum, then a maximum at a higher frequency. However, for 
shock-dominated clouds, the bremsstrahlung shows only a 
maximum at higher frequencies. 

Those are the most common types of continuum ob- 
tained from a single-cloud model. Very seldom the line and 
continuum spectra of the galax;ies can be reproduced by 
single-cloud models; on the other hand, multi-cloud models 
are generally required to explain the multiwavelengtli spec- 
trum and are obtained as a weighted average of the single- 
cloud models. The characteristics of some of the most com- 
mon single-clouds found in our previous analysis are listed 
in Table 1. Multi-cloud models lead to IR bumps wider than 
black body curves. Actually, in some compact objects, e.g. 
Mrk 3, Mrk 34, and Mrk 78 (Contini & Viegas 2000) dust 
emission is associated to a single velocity model. This is in 
contrast with the modelling of Circinus SED (Contini et al.), 
which requires clouds with different velocities contributing 
at different frequencies in the infrared range. That suggests 
that images of Circinus in the different IR wavelengths may 
show very different. 

Previous modelling of the representative galaxies with 
SUMA (e.g. NGC 5252, Circinus, NGC 4151) provided us 
with hints on which types of clouds are predominant in the 
NLR (Contini, Viegas & Prieto 1998a,b, Contini & Viegas 
2000). Adopting an initial grain size of 0.20 /um, it was found 
that Vs ranges between 100 and 1000 kms~^, no between 
100 and 1000 cm^^ , log(Fh ) between 8 and 13, and D 
between 0.01 and 10 pc. 

The IR bump results from summing up the SED con- 
tribution of dusty clouds with different velocities within the 
NLR, each of them peaking at the corresponding frequency. 
The maximum of dust emission occurs in the mid-IR for rel- 
atively high velocity clouds. The dust-to-gas ratio may vary 
from cloud to cloud. From the results obtained for Seyfert 
galaxies previously modelled by SUMA, we show in Fig. 2, 
the frequency range corresponding to dust emission associ- 
ated to a given shock velocity. A similar plot is shown in 
Fig. 3 for d/g. There is a tendency to find higher d/g for 
higher shock velocities. 



2.3 Sputtering of dust grains 

In models accounting for the presence of shocks, sputtering 
changes the initial grain-size distribution, creating a defi- 
ciency of small grains with radius agr < 50 A compared to 
their preshock abundances (Dwek, Foster, & Vankura 1996). 
The altered grain size distribution depends on shock veloc- 
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Figure 2. Vs adopted to model the IR data. NGC 7130 (dotted 
line), NGC 4151 (short dash), NGC 4051 (long dash), Circinus 
(short dash-dot), NGC 5252 (long dash-dot), NGC 4388 (solid). 
Mrk 3, Mrk 34, Mrk 78 (thin soUd) are selected from (Contini & 

Viegas 2000). 
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Figure 4. Effective temperature of dust as a function of the 
shock velocity for different agr (in /im). Dotted lines : agr=0.01 
//m, short-dashed ; 0.02 /im, long-dashed ; 0.1 /im, short-dash-dot 
: 0.2 \im., long-dash-dot : 0.5 \im., solid : 1 /im. Top panel : log 
Fh =11, bottom panel : log Fh =13. 



Figure 3. d/g ratios as a function of Vs adopted to model single 
galaxies. Symbols as in Fig. 2 



ity and on the density of the medium. The sputtering rate 

increases when the dust is in motion. At temperatures below 
10^ K, all the sputtering is caused by the relative gas-grain 
motion. 

When sputtering is present, the lifetime of grain survival 
depends on the distance downstream covered by the grains. 
In order to calculate the critical distance, we use the formula 
from ShuU (1978) which gives the rate of change of grain 
velocity behind the shock front: 

dVgr/dt=(vgr/2B)(dB/dt)-(7rOg^/mg,.)(pVg^)-Fpi/mgr 



Immediately behind the shock front the grain's gyrovelocity 
about the magnetic field B is Vg^ ~ 3/4 Vg. 

The first term on the right hand side is the rate of 
change of the grain's gyrovelocity calculated from the con- 
servation of the adiabatic invariant Vg^/B (Dwek 1981). For 
a "frozen in" magnetic field (see Contini & Shaviv 1982): 
(vg,/2B)(dB/dt) = (vg,/2n)(dn/dt), (B=Bo n/no) 
The second term contains the coUisional drag force resisting 
the grain motion (Dwek 1981) and the third term shows the 
plasma drag Fpj (ShuU 1978). Neglecting the plasma drag, 
if (vg,/2B)(dB/dt) < (tt a2,/mg,)(p v^,) 
the grains decelerate and sputtering is prevented. 

Considering that ds=dt Vgas is the differential distance 
corresponding to lifetime for evolution of the gas throughout 
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Table 1. The input parameters of single-cloud models 



model 


Ml 


M2 


M3 


M4 


M5 


M6 


M7 


Vs{kms-l) 


100 


500 


500 


500 


500 


700 


1000 


no (cm^'^) 


1000 


300 


300 


300 


300 


300 


1000 


log(Fh ) 




11 


11 


11 






11 


D (IQi^ cm) 


1 


1 


1 


1 


1 


1 


0.01 


SLgr ifim) 


0.2 


0.2 


0.2 


0.2 


0.2 


0.5 


1. 


d/g (10-")i 


1 


1 


10 


100 


10 


800 


10 



^ by number, corresponding to 4 10 * by mass for silicates 

a slab, and adopting a silicate density ~ 3 gr cm~^, we 

obtain : 8/3 (dn/n^) agr/mgas < ds, 

and integrating, the critical distance is given by: 

s (cm) > 10^^ &gr (pra) / a (10* cnr^) 

To perform the calculations a cloud is divided into a 
number of slabs (up to 300) where the physical conditions 
are considered homogeneous. The geometrical thickness of 
the slabs is automatically calculated by SUMA according to 
the gradient of the temperature downstream, in order to cal- 
culate as smoothly as possible the distribution of the phys- 
ical conditions throughout the cloud. The first slabs in the 
immediate postshock region show the maximum tempera- 
ture which depends on Vs. These slabs can be relatively large 
because recombination coefficients are lower the higher the 
temperature. The geometrical thickness of the slab closest 
to the shock front, S, is therefore determined by the model. 

Comparing 5 with the critical distance for models with 
low Vs(~ 100 kms^^), a preshock density n() = 100-300 cm~^ 
(corresponding to a density n ~ 1000 cm^'^ downstream), 
grains with agr=10~® cm are decelerated and hardly sput- 
tered if <5 > 10^'^ cm. For higher n (> lO" cm"^) and a 
stronger shock (Vs>300 kms~^), the critical distance is s 
~ lO"'"'' cm. Actually, the models show slab widths < 10"''' 
cm for Vs^ 700 kms~^. Since the physical conditions in the 
slab are assumed to be constant, sputtering may be efficient. 

For high velocities (Vs~ 1000 kms~^) temperatures are 
high. RecaU that da^r/dt oc (T/10^)^ ** (Draine & Salpeter 
1979), so small grains can be easily destroyed. Therefore, 
we adopt large grains and very small slabs (S < 10^^ cm) 
for grain survival. Narrow slabs improve the precision of the 
calculations. Such slab widths are also much smaller than 
the critical distances s which are ~ lO'^'' cm and ~ 10^^ cm 
for no=300 cm^'^ and no=1000 cm~^, respectively. 

The effective temperature of dust calculated for differ- 
ent shock velocities and different grain sizes is presented 
in Fig. 4. Notice that for Vs^ 700 kms^^only grains with 
agr> 0.1 /Ltm survive sputtering and contribute to heating 
the grains to relatively high temperatures. 



3 THE NEAR- AND MID- INFRARED 

CONTINUUM OF A SAMPLE OF GALAXIES 

The modelling of the near- and mid-infrared continuum of 
a selected sample of Seyfert galaxies presented by Alonso- 
Herrero et al. (2003) is presented in this section. The con- 
tinuum data provided by these authors has been subtracted 



from the light contribution of the background galaxy and are 
thus the best available data for analysis of the most central 
AGN emission. For modelling purposes, those objects pre- 
senting the widest wavelength coverage and minimum num- 
ber of upper limits were selected. For NGC 1068, the dataset 
from Rieke & Low (1975) is also included (open symbols). 
Objects which were previously modelled in detail axe omit- 
ted : NGC 4151 (Contini et al 2002), NGC 4388 (Ciroi et 
al 2003), NGC 5252 (Contini et al. 1998), and NGC 4051 
(Contini & Vicgas 1999). 

In this paper, we try to explain all the galaxies with pro- 
totype models. In all cases, clouds with Va= 100 km s"'^ are 
present and in most of them Vs= 500 kms~^clouds with 
different d/g also contribute to the fit. In a few cases, higher 
velocity clouds are also required. 

The results produced by multi-clouds models are shown 
in Figs. 5 and 6 for Seyfert 1.8-2 and Seyfert 1.5-1, respec- 
tively. The input parameters of the single-cloud models are 
listed in Table 1. For each galaxy, the single-cloud models 
used to fit the continuum are given in Table 2. The cover- 
ing factors of the clouds are related to the relative weights 
(see Table 2) which are used to obtain the best fitting multi- 
cloud model. Models Ml and M7 (Table 1) contribute sig- 
nificantly to the fitting of the SED in the near-IR and soft 
X-rays ranges, respectively. Their continua are shown in Fig. 
7. Model M3 is already presented in Fig. 1. 

To illustrate the results of multi-cloud models, in the 
panels corresponding to UM 148, NGC 3786, Mkn 335, and 
Mkn 530 , the single-model contributions are summed up for 
the dust emission and for the free-free emission separately 
and the best fitting SED is shown as a thick solid lino in 
Figs. 5 and 6. For the other galaxies the domains of the dif- 
ferent models arc well separated, so only single-cloud models 
appear in the figure. The summed SED would not change 
the picture. 

We note that in most cases, the data at the lowest fre- 
quency, 16 /im, should be taken as an upper limit, as it is 
largely contaminated by the host galaxy contribution due to 
the large ISO aperture. Therefore, the d/g for these clouds, 
which usually are associated with velocities of about 500 
km/s, is poorly constrained. In two galaxies, Mkn 334 and 
NGC 3362, the flux at optical frequency is also an upper 
limit. In NGC 7479 large d/g values (> 8 10"^^) are indi- 
cated by the data, in agreement with the upper limit in the 
optical domain. Data at higher IR frequencies are needed to 
confirm this d/g value. 

In brief, the infrared continuum of all the galaxies can 
be modelled by multi-cloud models obtained by the weighted 
sum of different single-cloud models. A fine tunc of these re- 
sults can be obtained by further modelling the emission-line 
spectra and the overal SED. Nevertheless, from the detailed 
analysis of a representative number of Seyfert galaxies (Con- 
tini et al. 1998a,b, 2002, 2003) no dramatic changes for the 
results of the IR-optical SED are expected. 

A comparison between Seyfert 1.8-2 and Seyfert 1.5 - 
1 modelling shows that IR-optical continuum fiuxes are gen- 
erally higher in Seyfert 1.5 - 1.0; moreover, the flux tends to 
decrease toward higher frequencies for some Seyfert 2 galax- 
ies. That behaviour is well followed by the models. For all 
types of galaxies, the continuum between 3 fim and 10 jj,m is 
mainly due to thermal bremsstrahlung by gas with relatively 
high density and heated by a low velocity shock (model Ml). 
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Table 2. The models in Figs. 5 and 6 



Galaxy Seyfert type models" relative weights'' CF'^ 



Mkn 334 


1.8 


Ml, 


M2, M3 


0, 


-4.1 , -4.6 


-8.9 


Mkn 573 


2. 


Ml 









-8.9 


TTM IzLfi 


1.9 


Ml 


M9 M*^ M4 

iVi^, IVi-O, iVi'i 


Q 


R -4 R -4 1 


-10. 




2.-1.8 


Ml 


M5 


Q 


-2.7 


-6.9 


IN VjV_/ OOUZr 


2_ 


Ml 


M4 


Q 




-10.5 


i> VT O / OU 


1.8 


Ml 


iVlOi i-Vl 1 


Q 


-4.5 -4.3 


-9 5 


vj \_j t y 


1.9 


Ml 


M2 


Q 


-3.75 


-9.1 


i'l vJ O OXJOO 


1.9-1.5 


Ml 


M3 M7 


Q 


-4.5 , -4.6 


-9.4 




2, 


Ml 


M3 


Q 


-5.1 


-9.5 


NGC 5273 


1.9-1.5 


Ml, 


M3 





-4.13 


-9.7 


NGC 5347 


2. 


Ml 









-9.0 


Mkn 471 


1.8 


Ml, 


M7 





-4.8 


-10. 


NGC 5674 


1.9 


Ml, 


M3, M7 





-4. , -5.3 


-9.5 


UGC 12138 


1.8 


Ml, 


M3, M7 





-4.4 , -2.1 


-9.4 


NGC 7479 


1.9 


Ml, 


M6 





-3.3 


-9.6 


NGC 7674 


2. 


Ml, 


M3, M5 





-5.3 , -3. 


-8.5 


Mkn 335 


1.0 


Ml, 


M2, M3 





-3.6 , -3.6 


-9.1 


Mkn 590 


1.0 


Ml, 


M3 





-4.6 


-8.8 


NGC 3227 


1.5 


Ml, 


M3 





-4.2 


-8.5 


Mkn 766 


1.5 


Ml, 


M3 





-3.9 


-8.7 


NGC 5548 


1.5 


Ml, 


M2, M3 





-4. , -4. 


-8.7 


Mkn 841 


1.5 


Ml, 


M3 





-4.1 


-9.1 


NGC 7469 


1.0 


Ml, 


M2, M7 





-3.95 , -4.6 


-8.4 


NGC 7603/ Mkn 530 


1.5/1.8 


Ml, 


M2, M3 





-3.8 , -4.5 


-8.9 



" The models are represented in Fig. 5 with the following symbols : Ml : long-dashed; M2 : solid ; M3 : short-dashed ; M4 : 
sliort-dash-dotted ; M5 : long-dash-dotted ; M6 and M7 : dotted. 

in logarithm 

The conversion factor (in logarithm) for each galaxy which accounts for the distance of the emitting clouds to the active center and 
for the distance of the galaxy to Earth. To compare the multi-cloud models with the data, the weights must be multiplied by the 
conversion factors. 

This indicates that dust emission in this region is relatively 
low, in contrast with most of the luminous IR galaxies where 
reemission by dust dominates at wavelengths even lower 
than 10 /im (Contini & Contini 2003). Clouds with Vs^ 
500 km s^^ contribute to the SED at higher frequencies. The 
relative weights of clouds with Vs> 500 km s~^ and clouds 
with Ve=100 km determine the slope of the multi-cloud 
continuum SED. In average, the ratios of relative weights 
for models Ml and M3 (see Table 2), which appear in al- 
most all the galaxies, arc higher for Seyfert 1.8 - 2 (3.7 lO**) 
than for Seyfert 1.5 - 1 (1.3 10**). These models correspond 
to velocities of 100 kms~^and 500 kms~^, respectively. Ml 
is shock-dominated while M3 is radiation-dominated. From 
the point of view of the unified model for the AGN, it is 
expected that high velocities clouds, reached by a more in- 
tense radiation, should play a more significant role in the 
spectra of broad-line objects, e.g. the NLSl galaxy Ark 564 
(Contini, Rodrigues-Ardila & Viegas 2003), since the lino of 
sight reaches deeper into the central region, where the high 
velocity clouds are. Accordingly, M3 appears in 9 out of 16 
narrow-line Seyfert but in 7 out of 8 broad-line objects. 

Interestingly, clouds having low Vs and high no, basi- 
cally model Ml, are present in all the galaxies. Whether 
these low velocity clouds are dusty or not, could only be 
found out by modelling their far-IR nuclear emission, which 
correspond to radiation by grains at rather low tempera^ 



tures. We notice that these types of clouds are difficult to 
trace in the nebular spectrum because they produce weak 
Hues. Moreover, they were not revealed by previous investi- 
gations of AGN spectra because complete datasets between 
1.1 and 10 /im were not available yet. 

The nature of these clouds is being discussed in detail 
in the case of the NLSl Mrk 766 - included in the Alonso- 
Herrero et al. sample- (Contini, Rodriguez-Ardila & Viegas 
2003, in preparation) . An analysis of the far-IR SED is there- 
fore needed to determine the d/g and thus better constrain 
the modelling. An exmaple of modelling the complete SED 
is given in the Appendix, where an analytical method to ob- 
tain a first representation of the far-IR to UV SED of Seyfert 
galax;ies based on our model approach is presented. For this, 
a more complete grid of models is given. 



4 CONCLUDING REMARKS 

In the last few years we have shown that composite models 
that account for the simultaneous effect of photoionization 
by a central source and by shocks produce a fair represen- 
tation of the SED of quite a large number of different AGN 
types, including Seyfert 1.5-2 galax;ies, NLSl, galaxies where 
the AGN coexists with staxbursts, active galaxies with very 
extended regions and luminous infrared galaxies. 
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Figure 5. The modelling of selected galaxies from Alonso-Herrero et al. sample: Seyfert 1.8-2 



In comparison with other models in the literature, our 
approach docs not resort to any particular dust/gas geome- 
try but explore the possibilities of two mechanisms expected 
to be present in any AGN, bremsstrhalung and dust repro- 
cessing, to reproduce the IR-optical SED. 

We discuss the different contributions to the IR contin- 
uum of AGN considering that both photoionization from the 
active center and shocks ionize and heat gas and dust within 
the surrounding AGN clouds. Collisional processes start to 
be important in the presence of shocks with velocities of 
about 100 kms^^. Therefore dust grains arc heated by colh- 
sion with gas to relatively high temperatures, depending on 
Vs. Moreover, the grains are sputtered crossing the shock 
front to different radius depending on the grain initial rar 



dius, the plasma density, and the slab width covered by the 
grains before being decelerated aud stopped by collisional 
drag. The geometrical thickness of the slabs is determined 
by the assumption that the physical conditions in the slab 
are constant, and are constrained by the choice of Vg, no, 
and agr suitable to the NLR conditions. 

Observations indicate that there is a gradient of densi- 
ties and velocities in the NLR, both increasing towards the 
active center. In the internal region of the NLR the radi- 
ation should be more intense. We found, however, that in 
most Seyfert type 2 objects the high velocity-density clouds 
are shock-dominated. If these clouds are closer to the centre, 
the effect of radiation must have been partly cancelled by the 
effect of an intervening dusty medium. Such dusty medium 





could be clouds with a very high d/g, which will emit in 
the far-infrared. On the other hand, the clouds contribut- 
ing to intermediate and type 1 Seyfert galaxies are mainly 
radiation-dominated. These clouds have a direct view of the 
central radiation as expected from the unified model. 

In this paper, the model approach is applied to a sample 
of galaxies of different Seyfert types for which the best esti- 
mates of the nuclear- stellar subtracted- optical-IR SED are 
available. With a limited range of parameter values - shock 
velocity, pre-shock density, ionizing flux, grain size and dust- 
to-gas ratio - which are dictated from our previous detailed 
modelling of representative Seyfert galaxies - the model re- 



sults arc encouraging as they produce a fair representation 
of the SED for all Seyfert type, in particular, they easily 
produce steep IR SED for broad-line type objects. 

Further, a first estimate of the dust-to-gas ratio for some 
of the clouds is provided. Still, far-IR data are needed to 
constrain better those values. Fine tunning of the modelling 
can be achieved by cross cheking it with the respective neb- 
ular spectrum and the SED at other wavelengths. Neverthe- 
less, the analysis of the optical-IR contirmum of a sample of 
galaxies provides a general view of the mechanisms produc- 
ing the emission, as well as on how they differ in the two 
types of Seyfert galaxies. 
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A comparison of Seyfert 1.8-2 with Seyfert 1.5-1 SED 
models confirms that continuum fluxes are generally higher 
for Seyfert 1.5-1 and that the flux decrease at higher fre- 
quencies can be reproduced for most Seyfert 1.8-2 galaxies. 
Recall that in multi-cloud models the single-cloud models 
are summed up adopting relative weights. We explain, there- 
fore, the flux decrease at higher frequencies by lower relative 
weights of clouds with higher velocities and/or reached by a 
strong flux from the active center. 

For all type of galaxies, the continuum between 1 and 
10 ^mis found to be produced by thermal bremsstrahlung 
emitted by gas with relatively high density, no (lOOOcm"''), 
heated by a low velocity shock, Vs (100 kms~^). These 
clouds produce weak lines and therefore their presence can 
mostly be accessed from the SED analysis. Because of their 
low velocity, we expect that dust re-emits in the far-infrared. 

To provide a first analysis of an AGN SED, an analyti- 
cal method is presented in the Appendix. In general, galaxy 
continua are explained by multi-cloud models, leading to 
extended IR bumps, which must be decomposed into single 
components. Indeed, the SED profile in the IR of a galaxy 
depends on the ratio between the dust reradiation peak and 
the bremsstrahlung emission in the IR, i.e. the dust-to-gas 
ratios in the single-cloud models, but depends also on the 
relative weights adopted to sum up the different contribu- 
tions, which are constrained by the line spectra. 

A linear fit of the bremsstrahlung emission (normalized 
to the value at 12 /xm) is adequate for \og{v) < 13.50, while 
the emission by dust can be represented by a black body 
(see Appendix). 
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Figure Al. The continuum of ARP 220. Symbols as in Fig. 7 



APPENDIX A: ANALYTICAL APPROACH 

To illustrate the procedure on how to fit the SED of a given 
galaxy using the present modelling, an example for the case 
of the SED of Arp 220 is given below. Arp 220 was chosen 
for two reasons: a rich dataset is available for this galaxy 
and a single-cloud model is enough to provide a good fit of 
the observed continuum from radio to ultraviolet. 

The observational data come from the NED 
(NASA/IPAC EXTRAGAL ACTIO DATABASE) (Benford 
1999, Carico et al. 1992, Chini et al. 1986, Condon et 
al, 1983, De Voucouleurs et al. 1991, Douglas et al. 1996, 
Dressel & Condon 1978, Dunne ot al. 2000, Ealcs ct al. 1989, 
Gregory & Condon 1991, Jarrett et al. 2003, Moshir et al. 
1990, Rigopolous et al. 1996, Soifer et al. 1989, Spinoglio et 
al. 1995, White & Becker 1992, Waldram et al. 1996, Zwicky 
& Herzog 1963) and from ISOPHOT (Spinoglio et al. 2002) 
The SED of Arp 220 (Fig. Al) shows a shape similar 
to that of model M7 (Table 1, Fig. 7). Bremsstrahlung 
explains the data in the radio, optical and soft-X ray after 
removing the background stellar contribution. In the IR, a 
single black-body type bump fits the NED (black squares) 
data as well as those from ISOPHOT (white squares). The 
SED can be accounted for with clouds that have Vs=1000 
kms"'^, no=1000 cm~^, and Agr = 1 /itm. The d/g ratio 
is determined by shifting the IR bump along the verical 
axis by about one order relative to the original model M7. 
Note that model M7 is originally calculated for d/g=10~^''*, 
leading to a too low IR bump (Fig. Al, thin dotted line). 
The IR data are explained, therefore, with d/g ~ 10~^^, 
which is relatively high, in agreement with the fact that 
Arp 220 is very dusty (Arp ot al. 2001) and is also classified 
as an ULIRG. 

In ordcn- to allow interested readers to use our model 
approach, a grid of single-cloud models is provided. In all 
models it is assumed that the shock acts on the edge oppo- 
site to that reached by the ionizing photons from the nuclear 
source, i.e., clouds are moving outwards. The following input 
parameters are kept constant as follows: the preshock den- 



sity, no = 300 cm^ , the geometrical thickness of the cloud, 
D = 10^^ cm, the shape of the ionizing radiation: a power 
law with power index equal to 1.5, the dust-to-gas ratio, 
d/g = lO"^"* by number, and the magnetic field. Bo =10"'' 
gauss. The elemental abundances correspond to the cosmic 
values. On the other hand, the shock velocity varies in the 
range 100 ^ Vs^ 900 kms~^, the initial radius of the dust 
grain in the range 0.01 < agr^ 1.0 //m, and the flux at the 
Lyman limit Fh is lO" and lO" cm"^ s"^ eV"^ (Table Al). 

The SED of galaxies, usually show the old stellar popu- 
lation contribution between about 3 lO^'* and 3 lO'^^ Hz. 
The data are nested between the black body curve cor- 
responding to a temperature of T= 3000 K (upper limit) 
and the bremsstrahlung (lower limit). Not always the lower 
limit is sharply defined. Nevertheless, the lower limit of 
the data in this range of frequencies roughly constrains the 
bremsstrahlung. It is, therefore, recommended to use a large 
range of frequencies, e.g., 10*-10'^ Hz to obtain consistent 
modelling in the difi'erent domains (radio, IR, optical, UV, 
and softX-ray). In Figs. A2 and A3 we present the SEDs 
calculated for models which appeax in Table Al. The model 
SEDs are calculated at the nebula. In order to compare them 
with the data the models must bo shifted by a factor which 
depends on the distance of the galaxy, the distance of the 
cloud from the center, and the covering factor (see Con- 
tini, Viegas & Campos 2003). Then the bremsstrahlung is 
adjusted considering the value at the maximum frequency, 
which depends on Vsand the maximum in the optical range 
which depends on Fh . The model is constrained by the IR 
peak frequency, which also depends on Vs. 

Notice that the models are calculated with d/g = 10~^* 
by number, so, the shift of the IR bump determins the 
d/g characteristic of the cloud. Complex shapes of the 
bremsstrahlung maximum and of the IR bump imply a com- 
posite model, so the other components must be chosen from 
Figs. A2 c A3 and shifted accordingly. 

In the following, based on these models, an analytical 
method is proposed as an alternative for using the SEDs 
given above. For each model, the free-free and dust emis- 
sions are normalized to their respective values at 12 /um. 
These values (Ntr and 'Ndust corresponding to the \og{h'Fv) 
at 12 yum) are listed in Table Al. Each of the models in 
table Al has a characteristics Fh , Vsand a,gr, for wich the 
corresponding black-body effective dust temperature Ta^st 
is given in Table A2. 

On the other hand, a linear fit of the logarithm of the 
normalized bremsstrahlung emission is obtained for 11.5 < 
log(i^) < log{i^iim)- In Table A3 the coefficients for the linear 
fit (A and B) and log(i^;im) are listed for each model. As seen 
in Fig. 1, for high frequencies, the behaviour of the free- free 
emission depends on the input parameters of the models and 
is in general characterized by two maxima and a minimum 
value. In Table A4, for each model the frequency and the 
i^Fy is given for those points identified as maxi, mini, and 
max2 (Fig. 1), as well as the fiux at log v = 15.0. 

Using the results listed in Tables Al, A2, and A3, 
the normalized gas emission can be obtained between 
11.5 < log(z/) < \og{i'iim)- The continuum from dust nor- 
malised (yd,y) is a black body with T=Tdust, i-c. yd,y = 
log(i/F^/i/i2F^j2) with log(i^i2) = 13.40, and Planck func- 
tion in photons cm~^ s~^: 
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Table Al. The logarithmic values of at 12 /xm for the 
bremsstrahlung and dust emission 
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Table A2. Effective dust temperature 
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Table A3. Linear fit to the logarithm of the normalized 
bremsstrahlung emission: y = A log(i^) + B 



Fh 


Vs 


log(i^Hm) 


A 


B 


lell 


100 


13.30 


0.77 


-10.42 


lell 


300 


13.10 


0.80 


-10.61 


lell 


500 


13.10 


0.80 


-10.60 


lell 


700 


13.30 


0.75 


-10.05 


lell 


900 


13.30 


0.76 


-10.15 


lel3 


100 


13.40 


0.82 


-10.98 


lel3 


300 


13.40 


0.76 


-10.15 


lel3 


500 


13.40 


0.76 


-10.14 


lel3 


700 


13.30 


0.81 


-10.76 


lel3 


900 


13.30 


0.77 


-10.24 



ydi,=31og(://i^i2)+log((exp(hz^i2/(kTd„st)-l) 
/(exp(hi//(fcTd„,t) -1)), 

where h/k = 4.80 lO^'^^. The thermal bremsstrahlung 
is obtained from the linear fit. 

Notice that for Vs= 100 kms~^the bremsstrahlung 
emission is completely independent of a^r. Nevertheless, the 
effective dust emission varies with a^rfor all the models. The 
effective dust temperature depends on how efficient is the 
sputtering. In all models the dust-to-gas ratio is kept con- 
stant to d/g — 10~^* . However, in the fitting procedure of 
a given galaxy, d/g shall be modified to account for the ob- 
served IR peaks in the mid- and far-IR SED, as described 
below. Note that the importance of dust relative to free-free 
emissions depends on the d/g parameter, particularly in the 
near-infrared range. 

As shown in the previous sections, the infrared-optical 
SED is usually due to two or more types of clouds, each 
one corresponding to characteristic velocity and grain size. 
In the following, the steps that should be followed to get a 
crude fit of the IR SED, are given: 

(a) the contirmum above log(:') ^ 13.50, is used to find 
the type of cloud (Vs and a.gr) explaining the observed data. 
The free-free emission dominates the continuum in that fre- 
quency range. Using Tables A3 and A4, the best model is 
chosen. 

(b) once the cloud characteristics are found, T^ust, 
defining the corresponding black-body curve, is obtained 
from Table A2. 
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Table A4. Characteristics of the bremsstrahlung emission for log(;/) > 13.50 



Vs 



maxi 



mini 



max2 



last 



lell 100 



all 



14.00/+0.033 14.28/-0.142 14.36/+0.344 15.00/-0.654 



lell 
lell 
lell 
lell 
lell 
lell 



300 
300 
300 
300 
300 
300 



0.01 
0.02 
0.10 
0.20 
0.50 
1.00 



13.90/+0.123 
13.90/+0.123 
13.90/+0.123 
14.00/+0.152 
14.00/+0.230 
14.15/+0.230 



14.20/+0.079 
14.20/+0.079 
14.20/+0.079 
14.25/+0.620 
14.37/+0.096 
14.30/+0.118 



14.36/+0.382 
14.36/+0.382 
14.36/+0.382 
14.36/+0.460 
14.48/+0.515 
14.36/+0.475 



15.00/-0.664 
15.00/-0.664 
15.00/-0.664 
15.00/-0.581 
15.00/-0.141 
15.00/-0.568 



lell 
lell 
lell 
lell 
lell 
lell 



500 
500 

500 
500 
500 
500 



0.01 
0.02 

0.10 
0.20 
0.50 
1.00 



13.70/+0.075 
13.70/+0.075 

13.70/+0.075 
14.00/+0.174 
14.00/+0.226 
14.00/+0.256 



14.00/-0.003 
14.00/-0.003 

14.00/-0.003 
14.30/+0.066 
14.38/+0.092 
14.38/+0.170 



14.14/+0.316 
14.14/+0.316 
14.14/+0.316 
14.36/+0.468 
14.48/+0.511 
14.48/+0.558 



15.00/-1.000 
15.00/-1.000 

15.00/-1.000 
15.00/-0.469 
15.00/-0.146 
15.00/+0.044 



lell 
lell 
lell 
lell 

lell 
lell 



700 
700 
700 
700 

700 
700 



0.01 
0.02 
0.10 
0.20 

0.50 
1.00 



14.10/+0.200 
14.10/+0.200 
14.10/+0.200 
14.00/+0.186 

14.00/+0.186 
14.00/+0.217 



14.20/+0.146 
14.20/+0.146 
14.20/+0.146 
14.30/+0.015 

14.30/+0.015 
14.34/+0.062 



14.30/+0.394 
14.30/+0.394 
14.30/+0.394 
14.36/+0.476 

14.36/+0.476 
14.48/+0.472 



15.00/-0.235 
15.00/-0.235 
15.00/-0.235 
15.00/-0.500 

15.00/-0.500 
15.00/-0.270 



lell 
lell 
lell 
lell 
lell 
lell 



900 
900 
900 
900 
900 
900 



0.01 
0.02 
0.10 
0.20 
0.50 
1.00 



14.00/+0.165 
14.00/+0.165 
13.76/+0.058 
14.99/+0.200 
14.99/+0.200 
14.99/+0.200 



14.20/+0.092 
14.20/+0.092 
14.00/-0.098 
14.32/+0.084 
14.32/+0.084 
14.32/+0.084 



14.32/+0.433 
14.32/+0.433 
14.16/+0.295 
14.40/+0.490 
14.40/+0.490 
14.40/+0.490 



15.00/-0.236 
15.00/-0.236 
15.00/-1.800 
15.00/-0.314 
15.00/-0.314 
15.00/-0.314 



lel3 100 



all 



14.36/+0.471 14.70/+0.412 14.85/+0.754 15.00/+0.698 



lel3 
lel3 
lel3 
lel3 
lel3 
lel3 



300 
300 
300 
300 
300 
300 



0.01 
0.02 
0.10 
0.20 
0.50 
1.00 



13.86/+0.128 
13.86/+0.128 
13.86/+0.128 
13.86/+0.128 
14.00/+0.235 
13.70/+0.020 



14.18/+0.058 
14.18/+0.058 
14.18/+0.058 
14.18/+0.058 
14.33/+0.204 
13.85/-0.023 



14.30/+0.400 
14.30/+0.400 
14.30/+0.400 
14.30/+0.400 
14.53/+0.512 
14.06/+0.410 



15.00/-0.504 
15.00/-0.504 
15.00/-0.504 
15.00/-0.504 
15.00/+0.066 
15.00/-0.602 



lcl3 
lel3 
lcl3 
lcl3 
lel3 
lel3 



500 
500 
500 
500 
500 
500 



0.01 
0.02 
0.10 
0.20 
0.50 
1.00 



13.76/+0.117 
13.76/+0.117 
14.36/+0.117 
13.76/+0.117 
13.76/+0.117 
14.00/+0.252 



14.00/+0.075 
14.00/+0.075 
14.54/+0.235 
14.00/+0.075 
14.00/+0.075 
14.38/+0.100 



14.33/+0.269 
14.33/+0.269 
14.60/+0.537 
14.38/+0.320 
14.38/+0.320 
14.48/+0.511 



15.00/-0.560 
15.00/-0.560 
15.00/+0.174 
15.00/-0.469 
15.00/-0.469 
15.00/-0.081 



lel3 
lel3 
lel3 

lcl3 
lel3 
lel3 



700 
700 
700 
700 
700 
700 



0.01 
0.02 
0.10 
0.20 
0.50 
1.00 



14.00/+0.174 
14.00/+0.174 
14.00/+0.174 
14.00/+0.191 
14.00/+0.191 
14.00/+0.221 



14.15/+0.135 
14.15/+0.135 
14.15/+0.135 
14.30/+0.040 
14.30/+0.040 
14.34/+0.075 



14.30/+0.412 
14.30/+0.412 
14.30/+0.412 
14.36/+0.459 
14.36/+0.459 
14.48/+0.468 



15.00/-0.310 
15.00/-0.310 
15.00/-0.310 
15.00/-0.409 
15.00/-0.409 
15.00/-0.228 



lel3 
lel3 
lel3 
lel3 
lel3 
lel3 



900 
900 
900 
900 
900 
900 



0.01 
0.02 
0.10 
0.20 
0.50 
1.00 



13.70/+0.056 
13.70/+0.056 
13.70/+0.056 
13.70/+0.056 
13.70/+0.100 
14.00/+0.182 



14.00/-0.508 
14.00/-0.508 
14.00/-0.508 
14.00/-0.508 
14.00/+0.044 
14.20/+0.122 



14.16/+0.282 
14.16/+0.282 
14.16/+0.282 
14.16/+0.282 
14.16/+0.350 
14.30/+0.394 



15.00/-1.400 
15.00/-1.400 
15.00/-1.400 
15.00/-1.400 
15.00/-1.000 
15.00/-0.322 
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Figure A2. Models. Symols as in Fig. 7 
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Figure A3. Models. 

(c) because the two curves are normalized by the re- 
spective values at 12 /im, the non-normalized curves must 
be used to verify the relative importance of dust and free- 
free, by shifting vertically the black-body component, which 
corresponds to finding the d/g for this particular object. No- 
tice that at this point the scale of the vertical Eixis is defined 
by the model, and only the shape of the curves are used. 

(d) The intermediate and far-infrared emission are then 
used to obtain other possible dust components (Tdust), 
which will define the types of cloud (Vg, agr) as well as their 
relative importance (d/g) by shifting the non-normalized 
curves vertically. This will also define the corresponding free- 
free components. 

(e) The final step is to determine the relative weights 
corresponding to each type of cloud to get the best fit to the 
observed flux values. 
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Symbols as in Fig. 7 



